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^ ■ Abstract 

(N 

We consider a confining Yang-Mills theory without Goldstone Bosons which could 

^ ■ describe the bosonic sector of the weak interactions. This model can be gauge invari- 

^ ' antly regularized on a lattice. A strong coupling analysis of the low lying bound state 

(^ . spectrum indicates that the vector isotriplet bound state (the right quantum number 

to represent the W-boson) could be the lightest state if the mass of the pseudoscalar 

isosinglet is raised sufficiently by the effect of the chiral anomaly (in analogy to the 



CK ' T)' mechanism of QCD). 



This work is a preliminary study in support to an intensive lattice Monte Carlo 
analysis of the model. 



1 Introduction 

Today, all the known fundamental interactions of elementary particles (strong, weak 
and electromagnetic) are described by the Standard Model (SM) [jll. The matter 
fields of the SM (quarks and leptons) and the Higgs boson are considered to be ele- 
mentary. They interact with each other by the exchange of gauge bosons which are 
also considered to be elementary. The structure of the theory is essentially deter- 
mined once the matter fields and their transformation laws under the local gauge 
transformations are specified. The SM has received a great deal of phenomenological 
support and the structure of the model has been confirmed to a high degree of accu- 
racy. 

In spite of the beautiful corroboration of the SM by experiments, we may naturally 
ask the questions: how elementary are quarks and leptons, the gauge fields and the 
Higgs boson? The possibility of discovering further substructure within the particles 
of the SM remains a viable option for the physics which lies beyond that model. 
The aim of introducing substructures is to construct a simple fundamental theory 
with a few degrees of freedom which should be able to reproduce the SM as an ef- 
fective theory. In this way one hopes to be able to solve three important theoretical 
problems of the SM: 1) the family replication of the matter fields, 2) the existence 
of too many parameters and 3) the fine tuning problem. Several models treat the 
quarks and lepton, the gauge bosons and the Higgs boson as composite systems. To- 
day, a conspicuous number of theorems exists which have ruled most of the existing 
models [0 and radically restricted the possibilities to construct realistic composite 
models. Two models have survived: the "Strongly coupled SM" (SCSM) ^ and the 
Yang- Mills theories without Goldstone bosons |Q. The aim of the SCSM is not to 
construct a simple fundamental theory beyond the SM, but to propose an alternative 
to the usual Higgs mechanism. However, the Yang-Mills theory without Goldstone 
bosons proposes to describe composite weak bosons starting from simple fundamental 
principles. 

To be precise this last model considers the photon to remain elementary and switched 
off. The weak gauge bosons W^ and Z then form a mass degenerate triplet. This 
model is a usual Yang-Mills theory, with gauge fields and fermion fields, but with 
a special choice of the degrees of freedom. The degrees of freedom are character- 
ized by the local gauge group, by the global isospin group and by the choice of the 
fermion fields (Dirac or Majorana spinors). To be plausible a composite model of the 
weak bosons has to reproduce the known weak boson spectrum: the lightest bound 
states have to be the W-bosons and heavier bound states have to lie in an experi- 
mentally unexplored energy range. The only possibility to have a Yang-Mills theory 
which reproduces the weak boson spectrum is to choose the degrees of freedom in a 
way that they naturally avoid bound states lighter than the vector isotriplet of the 
theory (which characterizes the W-boson triplet). This is possible if the unwanted 
light bound states which naturally show up as Goldstone bosons or pseudo Goldstone 
bosons in many models (like, for example, a pseudoscalar iso-nonsinglet, which would 



be the pion analogue of QCD) vanish by the Pauh principle, i.e. they are symmet- 
ric combinations of Grassmann variables. To obtain such a theory it is sufficient to 
choose the hypercolor gauge group to be SU{2), the isospin group to be SU{2) and 
the fermion fields to be Majorana spinors. 

In our work we will investigate the spectrum of this theory for two reasons: first, 
it is a very interesting project by itself to investigate the spectrum of a confining 
theory which exhibits a spectrum which is completely different from the QCD state 
spectrum and, second, this theory contains fundamental ingredients needed to de- 
scribe a realistic model of composite weak bosons. 

Because of the non-perturbative nature of the problem the model must be treated 
on a lattice in Euclidean space. We emphasize that the spectrum in Euclidean space 
is equal to the spectrum in Minkowski space. A lattice regularization a la Wilson 
is possible because the choice of the isospin group SU{2) allows us to replace the 
Dirac mass term and the Dirac-type Wilson term by a hypercolor gauge invariant 
Majorana expression. 

As a first step of our investigation we present a strong coupling expansion. We will 
explain in a formal way the technic that we have introduced to perform the strong 
coupling expansion. In particular we have developed an algorithmic way to easily 
identify and evaluate the diagrams of the expansion. This algorithm is very useful 
because the mass splitting between the lightest bound states show up only at higher 
order of the strong coupling expansion and to obtain a significant result one has to 
evaluate many hundreds of diagrams. 

From the result of this calculation we show that the vector isotriplet bound state of 
this theory is lighter than all other bound states provided that the chiral anomaly 
rises sufficiently the pseudoscalar isosinglet mass, in analogy to the rj' mechanism in 
QCD0. 

This work is a preliminary analysis in support to an intensive lattice Monte Carlo 
simulation where the mass of the pseudoscalar isotriplet is evaluated explicitely [0 . 
Outline of the paper: In section 2 we describe the model, in section 3 we explain 
the technic of our strong coupling expansion, in section 4 we classify the composite 
operators according to the CP eigenvalues and in section 5 we evaluate the spectrum 
of these CP eigenstates. 



2 The Confining Gauge Theory without Goldstone 
Bosons 

We start with a definition of our modelQ |^ . For pedagogical reasons, in tliis section, 
we formulate our theory using Weyl spinors in a second step constructing the Majo- 
rana spinors and rewriting the action in a QCD-like form. 



2.1 Definition of the fermion fields 

We consider a gauge theory whose fermion content is represented by a Weyl spinor 
F^^{x). Here a denotes the (undotted) spinor index (a = 1,2), A denotes the 
fundamental representation index of a global SU(2) isospin group {A = 1,2) and a 
denotes the fundamental representation index of the local SU(2) hypercolor gauge 
group (a = 1,2). Starting from this fermion field F the following hypercolor singlet 
operators for the low lying bound states can be formed: 

• Lorentz scalars: with fermion number +2 resp. -2 





S{x) = F{x)QF{x) 




S{x) = F\x)QF\x) 


Lorentz vectors: 


isosinglet resp. isotriplet 




V,{x) = F\x)a,F{x) 




y/(x) = F\x)a,T^F{x) 


Lorentz tensor: 






Bl,{x) = F{x)Qa,a,T'F{x) 




Bl{x) = F\x)a,~a.T'QF\x) 



(1) 



(2) 



(3) 

where {cr^} = {i, ctj) (fi = 0, . . . , 3) denotes the four component generalization of 
the Pauli matrices acting on the spin indices of the Weyl spinors in the Euclidean 
space-time and {T^} (/ = 1, 2, 3) are the Pauli matrices of the global isospin group. 
The matrix Q represents the antisymmetric matrix in spin, hypercolor and isospin 
space, which correspond to the Kronecker product of ia2,iT2,iT2 (the antisymmetric 
matrices in spin, hypercolor and isospin space, respectively) 

Q = ia2® iT2 ® iT2 



^In [Q an unfortunate sign error has entered the calculations due to a wrong factor i in the fourth 
and fifth term of eq. (3). To interpret the result correctly one has to interchange isosinglet with 
isotriplet. I thank H.Schleret for a discussion of this point. 



The notation F(x)QF(x), e.g., is a shorthand for the contraction of all the spin, 
isospin and hypercolor indices of the fermion fields with the matrix Q. Obviously 
among the composite fields (1), (2) and (3) there is only one candidate for a Goldstone 
boson, namely S, associated with a spontaneous breakdown of the fermion number 
transformation group F -^ e*"F by the condensate (S) = (S) ^ 0, which should arise 
in the confining phase. It is expected that the CP odd part of S (namely, S — S) 
acquires a mass by the chiral anomaly in analogy to the rj' in QCD. Notice that a 
Lorentz scalar isotriplet operator vanishes by the Pauli principle. As discussed above 
this last property is fundamental for any realistic electroweak composite model. 

2.2 The Wilson lattice action 

There is a SU(2) hypercolor matrix U(b) (in the fundamental representation) defined 
on each oriented lattice bond b. Our convention is that 

U{-b) = U\h) (4) 

An oriented path cj on the lattice is a set of bonds 

cj = 6i U ^2 U . . . U 6„ (5) 

such that the endpoint of hi is the startpoint of fej+i for 1 < i < n — 1. We can 
associate a SU(2) hypercolor matrix with uj by defining the path ordered product 

U{uo) = U{h)U{h)...U{hr,) (6) 

The cr— matrices are defined as follows: 

., ^ _ J a^ if b in +yu direction , , 

1 — a^ if b in —jj, direction 

We will denote by p C A an elementary plaquette of the lattice A and with dp the 
boundary of p. The complete Wilson fermion action on a lattice A written in the 
Weyl spinor notation is given by 

S = (3Y.TrU{dp) 

pCA 

-2k Y. {F\x)(y{h)U{h)F{y) + ^ [F{x)QU{b)F{y) + F\x)U{h)QF\y) 

The parameter P is related to the bare coupling constant by /9 = 4-; k is the hopping 
parameter and r (=1) multiplies the Majorana version of the Wilson term. The 
Weyl fermions are represented by the Grassmann variables F{x) and F^{x). Lattice 
points are denotes by x G A. (xy) C A represent the oriented bond between nearest 
neighbour pairs of lattice points. 



2.3 Compact notation 

Defining a^ = {—i, a^) one can easily check the following three relations 

ZT2Uib)Ur2 = -Uihf (9) 

iaaa^wa = a^^ (10) 

Qt/(6)^a^Q = U{hfa^^ (11) 

They can be used to rewrite the first part of the fermionic kinetic term in the action 
as a sum of two equivalent parts 

^ F\x)a{b)U{h)F{y) = 

b={xy) 

\ Y: F\x)a{b)U{b)F{y)+^- ^ F{x)Qa{b)U{b)QF\y) (12) 

b={xy) b={xy) 

This observation allows us to define the following four-spinor fermion field and to 
rewrite the action in a compact form, which is similar to the Wilson action of QCD: 

^-'^' ^ ( St') ) '^'' 

where the index a runs from 1 to 4 and indicates the two upper spinor components 
Ff^ and F^^ for a = 1,2 and the two lower spinor components F^^ and F^^ for 

a = 3,4. Its transpose is simply defined by tp'^ = (F-^, F^). 
In the notation (13) the isospin transformation F -^ e*"-^F is 



• < m 



T 
QTQ^ 



ip -^ exp <^ta \ ^ ^r^^T ( "^ 



The spin matrices are defined in terms of ordinary 7-matrices by 

Fi'M _ J r^ = r + 7^ if b in +/i direction , . 

1 r^ = r — 7^^ if b in —jjl direction 

The 7-matrices are hermitian 4x4 matrices, satisfying {7'^,7''} = 25^'^ which we 
choose as follows: 



and we also define 



75 = 7o7i7273 = f 01) ^^^^ 



With this choice of (/? and r(6) and with eq. (12) we can rewrite the action (8) in the 
following compact form 



S = l3Y,TrU{dp) 



pCA 



b=(xv) \ 



r(&)f/(&)(j QJ^iy) 



b={xy) 

-l:Y.^^i^)M^) (17) 



where the matrix A is defined by 

We note that the action (17) is invariant under the isospin transformation. In section 
4 we will go one step further and use fermion fields of the form 

^=(J °).a„dV' = .^(;? ;) (19) 

in terms of which the action closely resembles the Wilson action of QCD. 

S = {3Y.TrU{dp)-k J2 Hx)T{b)U{b)^{y) -Wij{x)ij{x) = 

pCA b={xy) X 

^ ScauAU) - ^^^i^ 

We emphasize that one has to be careful since, unlike in QCD, -0 and -0 are not 
independent fields. They are Majorana fields with hypercolor and isospin. 

3 Strong coupling expansion rules 

We will use functional techniques to give a simple derivation of Wilson's Feynman 
rules for the strong coupling expansion that we will formulate on the lattice. By 
coupling the various degrees of freedom to sources, we will obtain a compact formal 
expression for the Green's functions of the theory. Strong coupling expansions have 
been discussed in [7,8]. Here we analyse the changes needed for the Majorana case 
for an expansion requiring high order. 

3.1 Functional integration 

Quantization of the fermion degrees of freedom can be carried out by expressing the 
transition amplitude directly as a sum over all possible paths in field space connecting 



the initial and the final states. The partition sum is then 

[dip][dU]exp{-S{ip,U)} 



The sum over the paths for the boson degrees of freedom is a functional integral over 
the element of the gauge group U(b): 

[dU] = n dU{b) (20) 

hcA 

The sum over the paths for the fermion degrees of freedom is a functional integral 
over the anticommuting fermion fields (p{x): 

[d^] = n d^ix) (21) 

xeA 

Thus they are elements of a Grassmann algebra. Notice that in contrast to QCD 
here there are no independent conjugate variables ip. The integration (21) is defined 
as follows 

• Grassmann integration: In a n-dimensional Grassmann algebra, the n Grass- 
mann generators 9i, 62, ■ ■ ■ ,0n satisfy {9i, 6j} = (for i,j = 1,2,..., n). The 
symbols d9i, d92, . . . , d9n are defined by their algebraic properties {d9i, d9j} = 
{d9i, 9j} = and 

' d9i = 
d9i9j = 5ij 

From these rules one derives that for any totally antisymmetric matrix A 

f d9id92 . . . denexp{-{9, A9)) = VdetA (22) 

which is the essence of the Grassmann integration. We obtain also a formula 
with shifted variables 



d9id92 . . . d9nexp{-{9, A9) + {r], 9)) = VdetA exp( — [t], A'^r])) (23) 

where r] = [rji, ...,?7n) are Grassmann variables. 

The strong coupling expansion is most easily derived and compactly formulated 
through the introduction of external sources coupled to the degrees of freedom of 
the theory. In this approach. Green's functions are obtained by differentiating with 



respect to the various sources. Corresponding to the degrees of freedom t/"° (b) and 
'^&ai^) "^6 introduce the sources K"-"- (6) and A^^(a;), respectively^. 

Perturbation theory begins by breaking the action S into two parts: the unper- 
turbed action 5*0 and the perturbation action Sj. For our strong couphng expansion 
we choose 

So = -^-j:^^ix)Aipix) (24) 



Sj = pY.TrUidp)-k y: v^i^)(^ ?)r 

P b=(xv) V / 



ib)U{b)(l gV(?/) (25) 



A n-points Green's function can then be written as 

= ^J[dcp][dU]^il,^ix,)...^taA^n)exp{-Siip,U)} = (26) 

[dip] [dU]exp I -Soiip, [/) + E A^(x)y.(x) + ^ K{b)U{b) 



X,K=0 

where Z = J[dip][dU]exp{—S{Lp, U)}. The notation X^f, for example, is a shorthand 
for the contraction of the hypercolor, spin and isospin indices. The division by Z is 
equivalent to the omission of disconnected diagrams with at least one factor being a 
vacuum diagram (no external legs). 

We note that the unperturbed action depends on the fermion fields only. The func- 
tional integration over the fermion fields yields (using eq. (22-23)) 

[dip]exp I -So{<p) + Y^ X^{x)ip{x) \ = 
1 



[dip]exp\J2 



^if^{x)A^{x) + X^ixMx) 



= exp|-EA^(a;)^A(x)| (27) 

The functional integral over the gauge group gives us 

J[dU]explYK'''''{b)U'"''{b)\ =llDiK{b),K{-b)) (28) 



2 



A(a;) can be expressed by two sources associated to the two Weyl fermion fields F{x) and F''{x), 



^A _ ( Va,a{^) 



respectively: Ag^ ^ = _^'° , ( ) . Notice that the bar over rj simply denotes that rj and 77 are two 
different sources. 



where nb=(xj/) denotes a product over all nearest neighbour pairs of lattice points and 

D{K{h),K{-h)) = fdg exp{K""'(6)c/""' + ir""'(-&)(^"')""'} (29) 

Here dg denotes the Haar measure over the gauge group. Putting all this together, 
we obtain the expression for the n-point Green's function (we omit the hypercolor, 
isospin and spin indices for simplicity) 



X expl-J2^^{x)^X{x)\ X D{K{b),K{-b)) 



X,K=0 



' ' e-pl-S^(l,^]]xZ[X,K] 



6X{xi) SX{xn) { \5y 5K 
where the last factor Z denotes the generating functional 



(30) 

\,K=0 



Z[X,K] = - exp|-^A^(x)^-A(x)| x l[DiKib), K{-b)) (31) 

The expansion of the exponential of the perturbation part of the action exp{—Sj} in 
a power series generates the strong coupling series. The terms of this series may be 
represented by graphs similar to the Feynman graphs in ordinary perturbation theory. 
The graphical rules for calculating the Green's function can be read off from the equa- 
tion (30). In the next subsections we will present simple algorithmic Feynman rules for 
the two point Green's function of the form {ip{x)H(f{y)) or {ip{x) Miip{x)ip{y) M2(p{y)) 
where H, Mi and M2 denote some general matrices in spin, isospin and hypercolor 
space which produce hypercolor singlets and respect the covariance of the bilinear 
forms (fMjLp (for example Lp(f is not covariantQ). For the fermion string through the 
graphs along the path u we only need the familiar Dirac-fermion rules. But for the 
vertices at the external points x and y the rules are different from the Dirac-fermion 
ones, due to the Majorana nature of the fermion. 

3.2 Group integrals 

In equation (30) we are left with the integral D{K{b), K{—b)). This integral is a gen- 
erating functional for integrals of polynomials of group matrices trough the relation 



j dgg''^< . ..(,""< (^-i)^i^'i . . . [g-^f 
5 6 6 



6K''^<{b) ■ ■ ■ 6K^^<{b) 6K'''^i{-b) ' ' ' 6K^-^^'^{-b) 



D{K{b),K{-b)) 



(32) 

K=0 



^One can't contract dotted spinor indices with undotted ones. 

10 



The structure of the expansion is revealed considering the integrals over the group 
elements on the path u. A general discussion of integration over SU(N) group ele- 
ments is given by Creutz |^. For our purposes it is useful to know the following two 
results (together with the Haar measure normalization J dg = 1): 

dgg^^' = (33) 

dgg--\g-Y = ^S^h'Sa^b (34) 

We consider a closed path u which forms a loop from a point x G A to a point ?/ G A 
and returns to x 

uj = uj^ U uj~ 
uj^ : X I — > y 
uj^ : y I — > X (35) 

(the path can have any other closed form, it can be, for instance, the union of two 
disconnected closed paths: a loop from x to x and another loop from y to y). The 
graphs with non trivial group integral are given by the following algorithm: 

• The path: We draw the path u on the lattice A. We choose a direction of the 
path. Each bond 6 G w defines a link between two points. The direction of this 
link is defined by the direction of the path. 

• Plaquettes addition: In order to have a non vanishing group integral over a 
hypercolor singlet polynomial of group matrices we have to add plaquettes in 
the following way. 

A plaquette is bounded by four directed links. For each given path u we add 
plaquettes on the lattice A in such a way that there are 2n (n is some integer 
number) or zero links for each bond 6 G A and if there are 2n links at a bond 
b, half of the links are directed along +b and the other half along —b. 

• Vertices and bonds: We consider only the bonds with 2n links. This set of 
bonds is called the bond set and the respective set of start and end points of 
these bonds is called vertex set. The bond set forms a general graph. 

• Connected graph: We call graph G{uj) a connected general graph constructed 
from a path u with the previous three steps of this algorithm. The graph G{uj) 
is not unique, but the group integral of two different G{uj) of the same path u 
and with the same number of vertices is unique. 

• Group integral: We denote with Qab' a'b{G{uj)) the group integral of the path 
UJ for graph G{ll!). The hypercolor indices aa' and bb' indicate the hypercolor 
indices at the begin (a,b) and at the end (a',b') of the paths uj~^ and u!~ , 
respectively. 

11 



— Graph of type I: We call graphs of type I a graphQ constructed from a 
path u with the first three steps of this algorithm with n=l. All graphs 
^(a;) of type I of the same path u and with the same number of plaquettes 
have the same number V of vertices and B of bonds. The group integral 
for a graph of type I yields a very simple result: 

nab'a'b{G{uj)) = N^^~^-^'^5ay5a'b (36) 

A term N arises from contracting the hypercolor indices of eq. (34) at each 
vertex and a term 1/N arises from each bond with two links (see eq. (34)). 

— Graph of type II: We call graph of type 11 a graph which is not of type 
iQ. The group integral for a graph of type II has not a simple form and 
a general result like for the previous case. The group integral has to be 
evaluated separately following the methods explained in |^. 

3.3 Strong coupling Feynman diagram 

A strong coupling Feynman diagram D{G,uj) is formed by a path u and one of 
its graphs G{uj). The order of a term in the strong coupling expansion which is 
characterized by a diagram D{G,uj) is given by the number P of plaquettes of the 
graph G{uj) and by the length of the path |u;|. A diagram will give a contribution to 
the Green's function proportional to (— Z?)^ x fcl'^L The amplitude A{D{G,uj)) of the 
diagram is given by: 

A{D{G,u)) = {-l)'^nab'a'b{G{uj))^aa'w{uj) {-(3f&\ (37) 

where L is the number of internally closed loops formed by the path and S(a;) is the 
trace term arising from the trace over the matrices of spin and isospin. The trace 
over the hypercolor is obtained by contracting the indices aa' and bb'. This last trace 
term E is different from the analogous term that one expects from a strong coupling 
expansion with Dirac fermions and its evaluation requires care. Therefore we will 
discuss it in a separate subsection. 

3.4 The trace term 

In this subsection we consider only two-point Green's functions of the form: 

{^{x)M^{x)^{y)M^{y)) 
The normalized Green's function is 

{^{x)M^{x)^{y)Mi^{y)) = ^ J[d^][dUMx)M^{xMy)M^{y)exp{-S{^,U)} 



^SeeFig. la 
^SeeFig. lb 



12 



When we consider only the fermion integration, the normahzed fermion expectation 
value will be 

((^(x)Mv9(x)(^(y)Mv?(y))([/)^'^™" = ^ /[rfv^]v?(x)Mv9(x)y;(y)M^(y)ea:i9{-5((^, /7)} 

This can be written as a sum over paths forming fermion loops connecting x and y. 
There are two possibilities: 

• Isomultiplet: If the matrix M contains some isospin matrix T^ there is no 
contribution from separate loops at x and y, since T^ is traceless. Only loops 
connecting x with y are allowed. 

• Isosinglet: If the matrix M does not contain a T^ term the contributions from 
separate loops at x and y are allowed. If two loops are connected by plaquettes 
then they correspond to a chiral anomaly contributions. If the two loops are 
not connected by plaquettes (i.e. they do not form a connected diagram) then 
the truncated Green's function has to be defined as follows^ 

((^(x)Mv9(a;)(/7(2/)M(/?(2/))*™"""*"'^ = 

= {ip{x)M^{xMy)M^{y)) - {^{x)M^{x)){^{y)M^{y)) (38) 

For the moment we concentrate ourselves on the non-anomaly terms. We will discuss 
the anomaly later in section 5 . We have defined the closed loop path u> from a point 
X to a point y as the union of two distinct paths a;"*" and uj~ . Each of these two paths 
is the union of the \oo^\ bonds which compose them 

^~ = U (%•-!' %•) (39) 



X. 



where Xq = x, X|^+| = y- yo = y and y\^-\ 

We use the shorthand SX{x) = jA^ for the differentiation of the fermion sources. 
The contribution of the fermion along the path is expressed by the derivative of the 
fermionic sources at each lattice point on the path. The kinetic term, proportional 
to the hopping parameter k, yields terms of the following form 




Qaa'{u;+)Z = S\^{X) "^^ ^ r(cU+)x5,,, X ' ^ UA(|/) Z 




X{z)=0 forxy^Zy^y 



^Candidates here are the CP eigenstates 5'± — S±S. However, one can easily show that {S-) = 
(the contribution of a graph to {S- ) cancels with the contribution of the same graph with all the 
space direction reversed) which avoids spontaneous CP violation and (5*+) 7^ due to a spontaneous 
breaking of the U(l) fermion number transformation. The vacuum expectation value of the vector 
fields is zero, because a spontaneous breaking of the Lorcntz invariance should not occur. 
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.T,.,^ I Q \ ^,, ,_, .^ , . / 1 



ew{io-)z = 5X' (y) I ^ ^1 r(^-) X 6w X ( Q j '^-^^^^ ^ 



Hz)=0 forxytzjty 

(40) 



where 



Qaa'iuJ^) = Saa' Y[SX^{Xj^i)i ^ ^ j T{{Xj^i, Xj)) i ^ ^ j 6X{Xj] 






Q \t../ \n / 1 



Qwioo-) = S,y n'^A^(%-i) n 1 r((^.-i'2/.)) n n M^(%) (^1) 



j=i 



1 /-VNyj-i,yj/; 1 Q Q 



When we differentiate the exponential of the quadratic form X^A ^A in Z twice with 
respect to A we have 



A 



-1 



where 



^^a,a{^)^^^,b(^) X exp{-X {X)^-X{X)} = 

= {'^~AB,^p,a,+ ('^''^(^))AMA-'X{x))^^^^,) X exp{-X^{x)^X{x)} (42) 



AAB,a^,ab " ( Q Q j ^^^^ 



denotes the totally antisymmetric matrix in isospin, spin and hypercolor space. Eq. 
(40) is obtained using eq. (42) in (41)[]. The external source derivative ^j^ ,/ — j-yr 
associated to the gauge field can be viewed as an identity matrix 6cc' because of eq. 
(34) and the algorithm of sect. 3.2. The remaining terms 6aa' and 6bb' arise from 
contracting all the hypercolor indices of the 6cc' for each bond of the path. The 
remaining 6X{x) and 6X{y) in eq. (40) will contract their indices with the indices of 
the derivatives 6X{x) and SX{y) arising from the propagator {(p{x)MLp{x)ip{y)M(^{y)) 
when we substitute the fermion fields ip with the derivative 6X. Applying these 
derivatives to the quadratic form in Z we obtain 

^aa'bb'iuj) = Tr{6Xix)M6Xix) x Qaa'iuj^) X 6Xiy)M5X{y) x Qbb'iuj-)z}\^^^ (44) 
Using eq. (40) in (44) we obtain 

^aa'bb'i^) = ^aa'hv X 

xTr[5X{x)M5X{x)x5X{x)(^^ M r(a;+) (" J g ) '^A(?/)x 



x5X{y)M5X{y)x5X{y)(^^ \\v^u-)(\ M 5A(x) x Z 



(45) 



A=0 



^AU terms involving the second term in the bracket on the right hand side of eq. (42) will vanish 
because of the backtracking (one obtains gamma matrix multiplications of the form T{h)T{—h) = 0) 
or because they are proportional to a X{xj) which will be set to zero in eq. (40). 
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The result one obtains from eq. (45) may be very different from the analogous result 
for Dirac fermions. The difference is due to the fact that in the Dirac case there are 
two independent fermionic external sources (say A and A) while in our case there is 
only one. In the Dirac analogue to eq. (45) there are two dX and two 6X, and the 
contraction of the indices is evident. In our case we have four dX and also four times 
the same functional derivative which leads to a more tricky contraction of the indices. 
In eq. (45) one can contract the indices of the derivatives dX{x) and dX{y) in many 
different ways, because each 6X appears four times. The proper contraction of the A 
using an explicit notation leads to a sum of four terms 

xTr \r{u+)MT{uj-)M - T{uj+)AM'^ AT{uj-)M 

+T{uj-^)AM^ AV{uj-)AM^ A - T{uo^)MT{uj-)AM'^ A] (46) 



where 



--iiWill) («) 



4 CP eigenstates 

In weak interactions CP is a good quantum number but not C and P separately. 
Therefore, we perform a classification of the composite operators (1-3) according to 
the CP eigenvalues. The CP transformation of the Weyl spinor F is defined by 



^^^^-^ - - '-.F^: 



[x) = iT2i(y2r ' [Xp] 



(where xp = (t, —x)) which fixes the phase. In the compact notation we can rewrite 
the CP transformation in the following form: 

tp^'^^ix) = n^ixp) 

^^^(x) = ^ixphoT2 (48) 

We can build CP odd and CP even combinations^ of the composite fields (1), (2) and 
(3). These eigenstates can be expressed in a very compact form with the definitions 
(19). 

• Scalar CP eigenstates: The scalar CP even combination of scalars is defined 
by 

S+{x) = ^{x)iIj{x) (49) 

The pseudoscalar CP odd combination is defined by 

S'_(x) = i'0(x)75'?/'(x) (50) 



**Lorentz scalars S", vectors V^" and tensors T^'' are CP even if CPS = S, CP V = V^ and 
CPT^'•' = Tf,^ and are CP odd i{CPS = -S, CPV^' = -V^ and CPPt^" = -T^^, respectively. 
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• Vector CP eigenstates: The vector isosinglet CP even combination is de- 
fined by0 

V^{x) = '-'ilj{x)ri^^{x) (51) 

The isotriplet vector states is defined by 

V^'\x) = -^(a;)7^rV(x) (52) 

This state is CP even for 1=1,3 and CP odd for 1=2. Q 

• Tensor CP eigenstates: The tensors are defined as follows: 

B>"'^ = ij{x)a^"'T^ij{x) (53) 

and its dual 

Sf ^ = e^'^'^PBlp (54) 

where a'^^ = ^[j^, Y]- B^^^ has CP properties opposite to the vector isotriplet 
V^ ones. Its dual B'^^^ has CP properties opposite to B^^^ . 

The propagators of the CP eigenstates going from the point x to the point y along 
some path uj can be evaluated with the help of eq. (37) and choosing the appropriate 
representation of the matrix M in eq. (46). In the CP formalism we can write a 
propagator in the form 

{{v{x)M^{x))^v{y)M^{y)) = {{tlj{x)m{x))H{y)m{y)) (55) 

Inserting the matrices M relative to the CP eigenstates in the eq. (37) and (46) and 
using eq. (9-11) the result takes a simple form 

mx)M^l^{x))my)m{y)) = Y. J2 MD{G, u)) = 

= J2 Y. Cm X Tr{M^T{uj+)MV{uj-)) x (-l)^fi(G(a;))(-/?)^A;l"l (56) 

uj:x^y\^x G{ui) 

where Cm is some constant and Q{G{uj)) = ^ab'a'b{G{uj))6aa'Sbb' ■ Notice that for a 
graph of type I we obtain il{G) = A^(^~^) (see eq. (36)). 

^Notice that using eq. (10) we obtain 

V^' = F\x)a^'F{x) - F^ {x)Qd^'QF\x) = F\x)a^'F{x) + F\x)(J^'F{x) ^ 2FHx)at'F{x) 

-"^^Notice that the charged and neutral weak bosons are defined by 



w+ \ 


/ Wi+ lW2 


Ws - 


Wz 


w- 1 


V W^i - iWi 



The CP properties are: CPW-i{x) = W^(xp) and CPW^{x) = W-(xp). 
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5 Strong coupling expansion masses 

In order to compute the mass of the CP eigenstates we will consider the static prop- 
agators of these states. 



Glucit) = E((^((0,^))M^((0,f )))¥((*, ^))M^((t,:?))) 



(57) 



First we assume that xq = and yo = t and x = y (summed over all x) to interpret eq. 
(56) as a static propagator. The lowest order diagram representing a static propagator 
is a double fermion line in the time direction as shown in Fig. 2a. Since there is a 
non-zero probability for a transition from the static state to a more complicated 
state, the full propagator is given by a sequence of excited states connected by static 
propagators as shown in Fig. 2b. Each diagram D{G{u),u) of a general path cu 
can be viewed as a space-time process contributing to the excitation of the static 
propagation. These excitations renormalize the mass of the static propagator. 

5.1 The unrenormalized static propagator 

We denote with G'^^(t) a static propagator. For the CP eigenstates we will consider 
the operators given by the following set of M's.Q 



M 


flavor 


75 + 757o 


S^ 


^kT' + a^'T^ 


ykl 


757k 


yk 


(jkjrpi 


B^Ji 


1 


s+ 



We will use the shorthands 



5_, V^^, V\ B^^\ S^ 



to denote the flavor of the bound states and its corresponding matrix M . The non 
renormalized mass of a M state can be found from the property of the unrenormalized 
static propagator (for t -^ oo) 



G^{t) 



Cj^^ X exp^—rriQt) 



(5^ 



where C^j is a constant. The path ujq denotes the path of Fig. 2a. The diagram of 
this path contains no plaquettes. Computing the amplitude (37) of this diagram we 

"'^^ Because of the relations d^V^ = CxS- (when the anomaly is neglected) and df^B^"^ = C" x V^^ 
(where C and C" are some constants) and because of the Fourier transformation with zero momenta 
p in eq. (57) the static propagator mixes S- with V'^ and V''^ with B°'^^ . 
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have to evaluate the corresponding F matrices^ 

T{u+) = {vy = 2(*-i)ro 

T{uo) = (f 0)* = 2(*-i)f (59) 

The group integral of this diagram is fi(a;) = 2 and so we find from eq. (56) that 

Gf{t) = C^ X {Aky (60) 

The static unrenormalized propagator G^{t) is different from zero only for two CP 
eigenstates: 

M = S^, V^^ (61) 

From eq. (56) we find that for these bound states the non renormalized mass is 

mo = - log(4/c2) (62) 

The static propagator of the other CP eigenstates is zero. Their mass is mo = oo and 
remain very large even after renormalization. This result is a pathologic default of the 
strong coupling expansion. It does not describes any physics. The only information 
that one can obtain after renormalization is the spectrum of the states for which 
Go{t) 7^ 0. In the next subsections we will compare the masses of these two states. 

5.2 Renormalization of the static propagator 

Having computed the static propagator, we will now compute the renormalization 
produced by intermediate exited states. Fig. 3 shows us such an excitation, with 
initial point z and final point w, where the static unrenormalized propagator arrives 
and departs. For fixed z and w, we will sum over all intermediate exited states to 
obtain a total weight for the event. We denote this weight by D^\z,w). The full 
propagator takes the form 

G^'(x,2/) = E Gf{x,z,)D^'{z„w,)Gf{w^,Z2)D^'{z„W2)x 

Zi,Wi(i=l...n) 

x...xD'''{zn,Wn)Gfiwn,y) (63) 

where the points Zi and Wi are required to lie between x and y and to be ordered so 

that each Wi is "later" than Zi, which is itself "later" than Wi^i. This last eq. (63) 

represents the picture in Fig. 2b. 

In order to find the renormalized mass of the propagator, we consider the full static 

propagator 

G^{t)= J2 Gix,y) (64) 

yQ-xo=t 



^For the notation see eq. (14). We have fixed r=l which correspond to the case r=l in QCD. 
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Similarly we define^ 



and then (63) gives 



D^(t) = (Gf (t))-i E D'^x^y) (65) 

yo-^o=t 



Si,ti{i=l...n) 

~ e-'"«*ePM* (66) 

where tj and Sj are the corresponding time coordinates of Wi and Zj. So the renor- 
malized mass rrij^ is 

"^M = "^0 - Pm (67) 

The leading order oi pj^j can just be written in the form 

Pm= E {GS\t'))-'D'\0,w)= E {4er'D^{0,w) (68) 

u;=(t'>0,'u;) io={t'>0,-u;) 

To proof this last equation we have to insert it into the last term of (66) and to 
expand the exponential function containing pj^j, then we have to compare the result 
with the first term of (66). The excitation term D^^{z,w) is defined starting from 
eq. (37) in the following way: 

where the sum is over all closed paths u from z to w and return and over all graphs 
G{uj). A{D{G,u})) denotes the amplitude (37) of an intermediate excitation diagram. 

5.3 Results 

In order to calculate the renormalized mass of the two states V''^ and S^ to next 
leading order we have collected all the excitation diagram for all t' > and w in eq. 
(68) corresponding to the fc^ and f3^ perturbation order. Our results are valid in the 
asymptotic sense in a neighborhood of the origin in the {k, (3) plane. 

• Vector isotriplet mass: First we write the vector isotriplet mass rriyki up to 
fourth order in (3 and k. Our result is 

my.r = -log{4k')~3e + h'(3-h'p' + lk'f3'-^k'f3' 

2 4 8 16 

-IQk^ + lhk*(3 - ^k^(3^ + Uk^(3^ - ^k^/3' (70) 



"Notice that (Gj^(i))-i = ^e""* = ^(4fc2)-* 
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Pseudoscalar mass: We calculate the pseudoscalar mass neglecting the anomaly 
terms, since they arise only at order k^^ and are generally severely suppressed 
in strong coupling expansions |T^. These terms will be discussed later. In this 



calculation the pseudoscalar mass comes out to be less than the vector mass. 
This is not surprising in view of the suppression of the anomaly. In the next 
subsection the consideration of the anomaly terms will show that the pseu- 
doscalar bound state behaves like the rj' bound state in QCD and in particular 
that it acquires a positive contribution to its mass from the anomaly terms. 
Our result is 



153,. ^o 65,.^, 237 



36k^ + Uk^P - ^^k^(3^ + —k^P^ - =^k^(3^ (71) 

4 4 2 

5.4 The anomaly terms 

Here we will calculate the leading order anomaly contribution to the mass split- 
ting Atti^""™"" ^ = TTi*!* — m§_ where m*^^ denotes the pseudoscalar mass when the 
anomalies are turned on. We use the methods derived in section 3 to determine the 
contributions of the anomaly terms. The mass splitting is given by excitations which 
are characterized by two disconnected closed loops: 

^anomaly ^ fJ^ U C-J^ 

U'^ : X I — > X 
u^ : y\ — >y 

The diagram which corresponds to these two paths can be found by the algorithm of 
section 3.2. The amplitude of this diagram is given by the formula (37). The trace 
term can be evaluated by modifying eq. (45) to the new paths type. The result is 

y/ anomaly\ 

= Tr\5\{x)Mi5\{x) ^ 5\{x 




X 5\{y)M,5\{y) yi 5X{y) \ ^ ;' ]T{u;-){ ^ ;; ]6X{y)xZ\ (72) 

A=0 

For the pseudoscalar bound state S- the result takes the form 

S5_ (cj'^«°™«'y) = 4 [Tr (75r°r(cj^)) x Tr {^j5T^T (u^))] (73) 

To find all paths which yield a non trivial contribution to the trace term (75) one has 
to find all the loops a;* (i=l,2) for which 

Tr[^,T'T{uj')] ^ (74) 
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Clearly this last equation is true if T{uj^) contains a 75. This condition requires that 
uj^ covers bonds in any direction, since 75 = 7o7i7273- The reader can convince him- 
self that the shortest loops must have at least eight bonds. The allowed lowest order 
loops are plotted in Fig. 4. The spin trace (74) is the same for all these loops, up to 
a sign: Tr[j5T°T{uj')] = ±32. 

An excitation diagram which contributes to the anomaly mass splitting is composed 
by two of the loops in Fig. 4 (one for each end of the excitation) connected by pla- 
quettes to make a hypercolor singlet as required by the algorithm of section 3.2. The 
smallest excitations that can be made in this way are superpositions of two copies of 
the same loop, going in opposite direction. The lowest order of Am^ is given by the 
sum of all such diagrams. It is proportional to k^'^/S^ with a positive coefficient. 
The leading order contributions to Am^_ involving jS are proportional to k}-^(3'^ 
(n=l,...,4). These terms are generated by excitations made out of two of the loops 
shown in Fig.4, linked by n plaquettes. It is easy to check from these excitations that 
when (3 is turned on the effect is to increase the mass of the pseudoscalar bound state 
by 

V 2 8 32 / 

5.5 Conclusion 

We have obtained two classes of particles in this theory. 

• Finite mass: There are two bound states which have finite mass in the strong 
coupling expansion: The pseudoscalar bound state S- and the vector isotriplet 
bound state V''^ . The pseudoscalar bound state S- turns out to be lighter than 
the vector isotriplet bound state in the strong coupling expansion. It behaves 
like the rj' bound state in QCD and can acquire a larger mass by the effect of 
the chiral anomaly. 

• Heavy mass: The composite operators S+,B^^^ and V^ (see sect. 5.1) cor- 
respond in the strong coupling expansion to states with infinite mass. Such a 
behavior is known in QCD, where for example the axial vector meson gets an 



infinite mass in the strong coupling expansion [11|. In fact, in nature, axial 



vector mesons are much heavier than the vector mesons. 

The interesting feature of the use of Majorana fermions as basic constituents is that 
the spin one isotriplet bound states (having the correct assignment to be candidates 
for the weak bosons) can be the lightest bound state, provided that the chiral anomaly 
rises the mass of the pseudoscalar isosinglet sufficiently. 

We have used the results of this work as the starting point of a detailed investigation 
of the low lying spectrum of this theory by a lattice Monte Carlo simulation ^ . In 
particular we have numerically evaluated the contribution of the chiral anomaly to 
the pseudoscalar isosinglet mass. 
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Figure Caption 



1. a) A possible graph of type I. b) A possible graph of type II. 

2. a) A static unrenormalized propagator, b) An excitation of the static unrenor- 
malized propagator. 

3. Close up of an intermediate excitation. 

4. The four geometrically distinct loops which contribute to the trace (74). 
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